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Mean Potential "Variation at Sunrise and Sunset (Kew, 1912 and 1914). 

It is hoped that a detailed study of the synchronous variation of the 
atmospheric electrical quantities with the integral intensity of " strays " in 
wireless phenomena may be possible at Kew Observatory in the near future. 



On the Motion of Solids in a Liquid Possessing Vorticity. 

By J. Peoudman. 

(Communicated by Prof. H. Lamb, F.R.S. Eeceived March 28, 1916.) 

1. The following investigations were undertaken at the suggestion of 
Mr. G. I. Taylor as an extension of some of his own.* 

The equations of motion of any system of solids in a liquid which moves 
without vorticity were obtained by Lord Kelvin by a method based on that 
of ignorable co-ordinates in ordinary dynamics. They are of Lagrange's type, 
and were afterwards obtained by various writersf by an adaptation to hydro- 
dynamics of Lagrange's original method. 

Lamb's work leads to some hydrokinematical relations which we shall put 
in evidence below. 

For the rotational problems discussed in the present paper, an attempt to 
obtain a Lagrangian function has only been successful in one case. This 
case is that to which Taylor's principal results refer, viz., the two-dimensional 
motion, with uniform vorticity, which at infinity reduces to a rotation about 
an axis. That a Lagrangian function exists in this case is obvious from 
Taylor's results. 

^ An account of these investigations has not yet been published. The subsequent 
references are to a MS. which Mr. Taylor very kindly showed the author. 
+ See Lamb, ' Hydrodynamics,' 3rd Edit., Art. 136. 
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In other cases, instead of bringing into the equations the effects of the 
inertia of the solids, restriction is made to the finding of the generalised 
components of liquid pressure on the solids. These are denoted by Q^, 
where 

Q, = 2[[p^^S, (1) 

in the notation of the next section, jp denoting the pressure intensity. 

For two-dimensional problems with uniform vorticity, this is immediately 
effected by an extension of the ordinary pressure formula. 

For three-dimensional problems, constant uniform vorticity becomes impos- 
sible, and the motion of the liquid depends not only on the instantaneous, 
motion of the solids, but also on its past state. Throughout the consideration 
of these, problems, the motion is referred to axes which rotate with constant 
angular velocity. The determination of the pressure is first reduced to the 
solution of a Neumann's problem, and the method is then applied to- 
calculate the pressure components in any case of initial motion. The results 
are in agreement with the conclusions and experiments of Taylor. 

Two general problems of continued motion are then formally solved by 
means of infinite power series in the time. The first of these problems^ 
concerns the evolution of any given small relative motion inside a rotating 
envelope ; the second concerns the motion of a body without relative rotation 
in a liquid which extends to infinity and is there at relative rest. The- 
coefficients in the power series are derived from a certain sequence of 
functions of position, and each function may be found from the two imme- 
diately preceding ones by the solution of a Neumann's problem. 

The convergence of the power series is not here established, but for the 
special case of a spherical boundary the series have been shown to be highly 
convergent. The demonstration of this is reserved for a future paper, in 
which it is hoped to give a solution of the small motion of a sphere in a 
form which is numerically intelligible. The functions of position, above 
mentioned, are at present being evaluated in this case. 

To render the problems mathematically determinate in as simple a way 
as possible, motion at the outer boundary or infinity is prescribed. It is 
supposed that, if necessary, this prescribed state is maintained by external 
agency not here considered. 

Hyclrolcinematical Relations, 
2. Let X, y, z be Cartesian co-ordinates referred to fixed axes, and let 



(ih q.2, 



i 7«j 
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be independent co-ordinates of the Lagrangian type completelj specifying 
at time t the positions of the solids. Also, let 

be functions of x, y, z, t and ^i, ^2? . . . , £», such that (^r is the velocity potential 
of the irrotational motion in the case for which 

qr = 1, q^'s = {s^ t), 
there being no motion at infinity, if the liquid be not enclosed in an envelope. 
When the motion is irrotational, the velocity potential (j> is given by 

^ = i'l(j>l + q'2(l>2+'>'+q'n(l>m (2) 

but we shall also make use of the functions ^i, ... , <f)m when the motion is not 
irrotational. 

Let ^* denote the rate of increase of <jf) for fixed values of x, y, z, and q 
the speed of a particle of liquid. Then it easily follows from the work of 
Lamb (loe. cit.) that 



\i^* 



■-lo^^^dS^j^t 



^^ ' dn 







9'^'-,;S-i^2|U|^^S, 



m 



dq, 



(3) 



where djon denotes differentiation along the outward-drawn normal to the 
surface of a solid, where the surface integrals are taken over the surfaces of 
the solids, and the summation extends to all the solids. 

To simplify matters it is here supposed that the liquid extends to infinity 
and is there at rest. 

Now on transforming the surface integrals on the right-hand side of (3) 
into volume integrals taken throughout the liquid, and then differentiating and 
re-transforming to surface integrals, we see that 



(a/ ^ 



h 



dn 



dS 



J L Sn ^ ^ dfi \dx dx dy 3y 9^ dz j dnj ' 

(4) 



1 



^ x\U^^dii = %\\ 



IrS 



d(f), 



dqr^n ^^ dn 



dS, 



(5) 



so that (3) gives 

2L^ 



A-— ./^^ ?^^ ^^^ii] ^^\^i 7Qi V 

r In.. ^-+9^ dy dz dz/'^dqrJdn 



\ox ox 



q^^4^dS = 0. 



jCigam, 



dn an an 

qr dqr oqr ' 



(6) 



dqi 



dqn 
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3^ d^r I 3^ 3^r , d(f> d<f>r __ / 3^1 3<|>r 3^1 3^r ,_ 9<^1 9^r \ • , 

a^ 0^ 01/ cy az az \ ox ax oy ay oz az J 

, [d(f>n d(f>r , 3^^ d^r , 3<f>n 3<^ A ^ 

\ 3^ 3^ dy dy dz dz J 

\ax ax -^ / \dy oy I \cz ^ dz 

m that, on equating to zero the coefficient of q's^t in (6), we ohtain 

'3<^^ __ dis\ _ l d<j>r 3^g d(l>r B(f)s , d(f>r 3<^g \"\ 3^ 

dqs dqrl \dx dx dy dy dz dz J J dn 



2 

n I > 






qt Sqrl \ dx dx dy dy dz dz JJ dn 

dS = 0. (7) 



+2 f^^+^-if^+^^\ ^ 

L 3^ 3^ dy dy dz dz J dn 



This relation remains valid if either s or i is equal to r or if they are equal 
to one another ; in both these cases it takes the same simpler form. 

As an illustration of the relation (7)> let us take the case of a single solid 
with qi, q2, q^ the Cartesian co-ordinates of a point fixed in it. Then on the 
surface of the solid we have 

_ — — ^^ — _ — ^^ __ if(,^ — —-- __ /f,j 
an on on 

ij m, n denoting the direction cosines of the outward-drawn normal, and 

-A-— — ^ 3 3 3 __ 3 

cqi ox aq2 ay oqs az 

when applied to <^i, ^2,^3. Taking r = 1, s = ^ = 2, in (7) we have 



f 



A dx dy J \ Sc dx dy dy dz dz 



the truth of which is easily verified by Green's theorem. It is interesting to 
notice the connection between this relation and those used by Prof. Lamb* 
to derive from the pressure formula the equations of motion for this 
particular case. 

TwO'JDimensional Frohlems with Uniform Vorticity. 
3. Let ^ denote the stream function and m the vorticity, so that 

^i^f = 2m, (9) 

•^ " On the Forces Experienced by a Solid Moving through a Liquid," * Quart. Jonrn. 
Math.,' vol. 19, p. 66 (1883). 
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whore Vi^ denotes v^ldx^-^-d^jdi/, x, y being Cartesian co-ordinates in the 
plane of the motion. 

Since ViV = ^^^ (10) 

we can find a function (^*, such that c^', ^/r* are conjugate functions, i.e. (f)' +iy}r' 
is a function oi x-\- iy. 

Now the equations of motion may be written in the form 



3? 6 



— 2oiV = 



31^ 

vx 



^ 



where P = p/p + ^(f + H, 

n being the potential of the external forces. 

Since 



Y 



(ii> 



(12) 



OX 



we have 






ct 



w 



dt dx dy 

and the equations of motion may be written 

— ^ — ^^ "5^ = ~ ^ ' 

ox ox ox 

_8^_2 ^= -^ 

"" 8;// 8?/ ' 



A 



Y 



(13) 



By 



V 



From these we immediately obtain the integral 

^//) = (/);+2a).;r— ir/-0, (14) 

it being supposed that the arbitrary additive function of t in (j)' is chosen so 
that (14) is valid as it stands. 

From (14) the particular results of Taylor, for a circular cylinder in a liquid 
whose motion at infinity is one of shearing, may be very easily reproduced. 

Now let us refer the motion to axes which rotate with constant angular 
velocity equal to ay. Let ^i\ v refer to the motion relative to these axes, and 
let y\r' be the relative stream function, so that 



dy 



^ t / 



v'='^, vi Y' = 0. 

OX 



1(15) 



The equations (11) remain valid with accented variables, if we take, instead 
of (12), 

V =p/p + h'' + ^-W'''\ (16) 
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q' denoting the relative speed of a liquid particle and r the distance from the 
axis of rotation. Taking <^' so that <^\ i/r' are conjugate functions, we obtain, 
as before, the integral 

In this case ^' = (jbi^'i + . . . + *^n9,'n' 

The foregoing results of the present section were obtained from the lectures 
of Mr. Herman, who quotes a paper by Prof. Love* as virtually containing 
tliem. 

4. Taking the case of rotating axes with relative rest at infinity, let us obtain 
the pressure components on unit length of a system of solid cylinders. "We 
have, from (1) and (17), 

071 



where ^^ = S 

P 
If 



Mr 



p 



= 2a,2[t'^is, (19) 

'^=Wt\r-^^Js. (20) 

5. We will now obtain the equations of motion for the problem just 
considered when O = 0, by the methods of Prof. Lamb. 

We consider only the dynamics of unit length of the system parallel to the 
axis of rotation. 

Let-^, 7] be the co-ordinates of a filament of mass m, whether of liquid or of 
one of the solids, and let X, Y be the components of force acting on it. 

We have 

:Sm(f A|^+^**A7;) + 2toSm(f A7;-7;*A^)-ft)2Sm(|A^4-^A7;) 

= :S(XA|+YA7;), (21) 

where Af, At; are the components of a " variation," and the summation 
extends to all filaments of the system. 

Now, just as in the case considered by Lamb, 

Sma AH, A,)=|^^(g-)-.g--|A,.+ ...+ |-^^(^^^ 

(22) 
where T' is the total kinetic energy of the relative motion, and 

t (XA|+ YAt;) = QiA^i + . . . + Q,A^„ (23) 

■^ " On the Stability of Certain Vortex Motions," ' Lond. Math. Soc. Proc.,' vol. 25, 
p. 18 (1893). 
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where Qi, Q2, ..., Qw, here denote the generalised components of the external 
forces applied to the solids only. 

In the present problem we have also 

t m (^A^+vAv) = i (~ A^i + . . . + |t A^,) , (24) 

where I is the variable part of the total moment of inertia of the system 
about the axis. It may be regarded as formed from solids occupying always 
the positions of the actual solids but of density G-—p, where a is the actual 
density of the solids, in general, of course, variable. 
Further, the contribution to 

2Sm(fAry-77-Aa (25) 

from the solids only, is 

where Mo' is the moment of momentum, about the axis, of the relative motion 
of the solids. The contribution from the liquid is 



-2p{ 



'l^'+t^f)''^. 



taken over all the liquid, which is equal to 

2p 2 1 '^An els, 

in our customary notation, where A?^ denotes the component of variation 
along the outward-drawn normal to the boundary of a solid. 

Now we might express the motion of the solids by means of a stream 
function ; and although this would not be a harmonic function, it could be 
taken equal to y\r at the surfaces of the solids. If, then, the density of the 
solids were uniform, and equal to po, their contribution (26) would be 

— 2po 2 1 -^ A'?^ ds. 

From this we see that the actual contribution from the liquid is equal 
and opposite to that of a set of solids of uniform density p, occupying always 
the positions of the actual solids. 

The total contribution to (25) is therefore 

where M' is the moment of momentum, about the axis, of the relative 
motion of solids, occupying always the positions of the actual solids, but of 
density cr—p. 
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The equations of motion may therefore be written 

at \oq ff clqr 

where E = T'4-«M' + -|-«2I, (28) 

the functions T', M', I, Q^ being formed as described above. 

Using this to obtain the generalised components of liquid pressure on the 
solids, we have 

«-- {l(|f)-|f} ■ (^«> 

qr''=-Wh> (31) 

where T/ is the kinetic energy of the relative motion of the liquid, Mi' and 
Ii are the moment of momentum of relative motion and moment of inertia 
respectively, of solids occupying always the positions of the actual solids, but 
of density />. 

The equivalence of (20) and (31) is easily seen, while the equivalence of 
(18), with 12 = 0, and (29) follows from the relations of § 2. To show the 
equivalence of (19) and (30), let us express the possible motions of the 
solids by means of stream functions i|ri', '>^^\ . . ., t^^', so that on their boundaries 
we have 

9s dn ' 

III/ \ Vl£ f J V(£i' \ O^j' 0(lr/ 



^ ]]\oy dx dx otj J 



the summation and integration extending over the solids only, and this is 
equivalent to 



2pX 



ylr --^ as 

^ OS 



= 2,SJt'f.. 



Three^Dimensional Problems^ 

6. We shall use Cartesian axes Qxyz rotating about O2; with constant 
angular velocity <», Let % v, w denote the components of the velocity at 
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.^, y, ^, relative to these rotating axes. Then the hydrodynarnical equations 
may be written in the form 

8P 



die 
dv 



2(v^—W7])—2cov 



dx ' 

3P 



>■ 



no 



dt 



2{%(/rj—v^) 



ap 



dz ' 



^ 



(32) 



where 






9 



V 



dib du 



X-, 2^ — 7-- 



9'^^ 3 



u 



(33) 



9^ 9^' ' dx 9y' 

and P = pIp + |^H O - ^^^ (^2 ^ ^2)^ (34^ 

q now denoting the relative speed. 

For the generalised components of liquid pressure on the solids, we have 



P 



= t[[{V-^lif^^^lc,\x^^if)}^^d^, 



or 



^r 



dcj^r 



9li 






(35) 



where Ii denotes the moment of inertia, about the axis of rotation, of solids 
occupying always the positions of the actual solids, but of density p. 

We shall now show that, for a given state of motion, the determination of 
the function P may be reduced to a Neumann's problem. 

On using the equation of continuity 



9?^ 3'y dvj 
dx dy dz 



0, 



we easily obtain from (32) and (33) at points in the liquid 

where v^ = 8^/9^2 + '^jdy'^ + 9^/9^^. Also, at a boundary of the liquid, we 
have, from (32), 

9P idu , dv , dia 

— 2 {l{v^--' W7)) + m (w^—itO + ^'- C^?—'?^? )} — 2co (lv—07mi\ (37) 

The whole of the right-hand sides of (36) and (37) is supposed known, 

except 

idu , dv , dvj y^^. 

iTTT + m^ + n-^, (38) 

at at dt ^ -^ 
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We shall now show that this also is determined by the instantaneous state 
of liquid velocities and the velocities and accelerations of the boundary. 

Let F {x, y, z, ^) = 0, 

be the equation of a bounding surface, which may be the surface of one of the 
solids. Then, if 

D _ a , a , a , a 

D^ ~^ dt doc dy dz 
we shall have 

^l\x,ij,z,t) = Q for Y{x,y,z,t) = 0, (39) 

where s is any integer. On taking s = 2 we have 

'D /oF , aF , aF , aF\ . ,,^, 

Dt \ at ox ay dzj 

on the surface. Expanding this, we see that the time derivatives of u, %\ w 
only occur in the form 

a^ 'bx dt dy dt dz ' 

and thus have (38) expressed in the form stated. 

The right-hand sides of (36) and (37) being known, the determination of 
P is reduced to a ISTeumann's problem. 

7. We shall now apply the methods of the preceding section to determine 
the pressure immediately after the sudden generation from relative rest of any 
possible motion of any magnitude. 

Let (j) denote the auxiliary function defined as in § 2, which would always 
be the velocity-potential if the axes were at rest and the motion irrotational. 

Now the initial relative motion will be irrotational, and (p will be its 
velocity-potential. Thus we have, from (36), 

V'i* = 0. (41) 

while it follows, from what was said in the preceding section, that 

-^^=-^4L + 2co(l^^-m^^). (42) 

on on \ oy oxj 

We therefore have, initially, 

P=^- + 2a)P^ (43) 

where P' is determined by 

in the liquid, and 



V'P' = 0, (44) 

.^^ = l^A^,n% (45) 

dn dy dx 



at the boundaries. 
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For the pressure components we have then 



where 



P 



t\\W-lcf-^fPd^ 



= 2co 2 



J 



2a> S 



— 2a) S 



8P' 



dn 



(j)r(l 



J J 



9^ 



0x1 



^r 



/// ':^ 



-— __1...2 



oil 



(47) 



(48) 



Here Qr is what the whole corresponding pressure component would be if 
the axes were at rest, while Q^''' is, as in two-dimensional cases, equal to 
the corresponding centripetal component of mass-acceleration of solids 
occupying always the positions of the actual solids but of density p. We 
may write Q/' in the form 



where 



2cop 

Hr,s 



/3r,l^*l + /3r,2^*2+ > • . + /5r,7i2*W' 



(49) 



-S 



4>r \l 



d6s 9<jf) 

dij 



ox 



dS 



-- S N)r sin 6 



ds 



dS. 



(50) 



where d/ds denotes differentiation along an arc of the boundary which lies 
in a plane perpendicular to the axis of rotation, and the angle between the 
outward-drawn normal and this axis. We notice that ^r,s = — ^s,r' 

The coefficients ^^.,5 are all that have values peculiar to the present problem. 
It is interesting to compare them with the corresponding coefficients in the 
gyrostatic components of the mass-acceleration of solids occupying always 
the positions of the actual solids, but of uniform density. These coefficients 
are given by 



Pr,s — 2 



dm, dm, 



dq'rdqs Sq'sdq, 

where Mo is the moment of momentum about the axis of the relative motion 
of solids occupyiDg always the positions of the actual solids, but of unit 
density. 

8. As an example we will calculate the coefficients ^r,s for the initial 
motion of an ellipsoid of revolution, with its axis parallel to that of rotation. 
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Let the co-ordinates of its centre be qi, q^, q^ ; then we require the functions 
4>ii <i>2y <^3 on the surface only. 

Let us first consider the case of the sphere of radius a. Take spherical 
polar co-ordinates r, 6, ^y with the pole at the centre of the sphere, the axis 
^ = parallel to Oz, and the plane ;)^ = parallel to that of y = 0. Then 
we shall have 

^1 = Aa sin 9 cos ^, <^2 = Aa sin d sin x^ ^3 

where A = |, while sin 9 d/ds = djacl^^ Therefore, 
A,3 = -A%3 



Aa cos 6, 



(51) 



(52) 



sin^ cos^ ^ cl6 cl')(^ = — fTra^A^, 

while ^1,3 = ^2,3 = 0. 

Let us next consider a prolate ellipsoid, of polar and equatorial radii 
a and h respectively. Take co-ordinates |, t], ;)^, such that if x^^, i/o> ^0 be the 
Cartesian co-ordinates of the centre of the ellipsoid, 

x^Xq = c sinh ^ sin tj cos x> P—'l/o = ^ sinh ^ sin rj sin x^ 

z — Zo = cosh f cos rj, 
where c^ = a^—¥. 

On the surface of the ellipsoid we then have* 

(^1 = Ah sin 7] COB ^^ ^2 = Ah sin t] sin %, <^3 = C& cos t], (53 



where 



A= |log{(l + .)/(l~.)}^./(l--e^) 

ilog{(l-i-6)/(l-e)}-6(l-26^j/(l-62y 



(34) 



^ being the eccentricity of a meridian section, so that ¥ = a^(l— e^), and is 
a similar constant, but is not required in what follows. Therefore 



/3l,2 = 






(f)i sin 6 -^ d8 



r2TT 



z:= -a Jx = Q ■^ 



A^a&2 



tr 



'27r 



sin^ 7] cos^ % r??7 d^ 



— 4.7rah^A^, 



(55) 



while ^1,3 = ^2,3 = 0- 

Considering an oblate ellipsoid of polar and equatorial radii h and a 
respectively, similar analysis shows that 

^1,2 = -^7ra%A^ A,3 = /32,8 = 0, (56) 

___e(l — e^)—- (1 — 6^)1 sin"i 6 
(3 (1 + ^2) — (1 — - e^)i sin" ^ e * 

* See, for instance, Lamb, 'Hydrodynamics,' 3rd Edit., p. 135. 



where now 



jfx 



(57) 
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"N'ow let us compare (52), (55), and (56) with the corresponding coefficients 
j3rj in the mass-acceleration of an ellipsoid of uniform density. We see that 
the ratio ^1,2! ^1,2 is in each case A^. IsTow in (51) A = ^, in (54) 

A — ^1 as e — ^ 1, 
while in (57) 

A — >- \ as e — ^ 0, 

A — ^ as 6 — ^ 1. 

Thus for the sphere the ratio ^i,2//3i,2 is ^ ; as the prolate ellipsoid tends 
to take a spindle shape the ratio tends to unity, which is its value for two- 
dimensional motion, while as the ohlate ellipsoid tends to take a disc shape 
the ratio tends to zero. These results completely bear out the general 
reasoning of Taylor. 

9. When the relative motion is either irrotational or small, the equations 
(32) become 



"N 



y 



du o 9l^ 

Ot OX 

■|;+2««=-^, y (58) 

a small motion being one in which squares and products of velocity - 
components may be neglected. It immediately follows that the only such 
motion which can remain '' steady '' is one in which 9P/3>: = 0— ^.e., a two- 
dimensional motion perpendicular to the axis of rotation. 

As a first problem in the discussion of motion over an interval of time, let 
us consider the motion of a liquid inside an envelope which is fixed relatively 
to the rotating axes, on the assumption that the motion always remains 
small. 

We have the equations (58) with 

P rr: -p/p + 0-~|6e)2(<:rH2/'), (59) 

and, at the boundary, 

lu + Qiiv-hniu = 0. (60) 

Suppose the following expansions to be permissible, at least for a certain 
range of ^^ : 

P==2a)ip,^P-^. (62) 

6- == "5 1 
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Here tig^ Vg, lOg and P^ are functions of x, y, z only, and %i^, %, %d^, which 
give the initial state of motion, are supposed known. 
On substituting in (58) and (60) we obtain 

ap. 



"A 



Ih^l — 'Os 



'^5+1 + % = 



^^.5+1 



'dx 

ap, 

8P, 

'hz 



y 



(63) 



^ 



with his'i-mvs'i-mvs = 0, 

at the boundary, for 5 = 0, 1, 2, ... . From these we derive 



(64) 



OX dy 
m _L ^, _ SP^ , oPs-i 



■^ 



dx 



Ws+l-^lVs-i = 



cly 

8P. 8P.-2 



y 



dz 



a 



''Z 



(65) 



^ 



for 5 = 2, 3, 4, ..., while these equations remain valid for s = 1 if we replace 
dV-ijdz by —Wq. 

From (65) we have then, inside the liquid, 



V^Ps = — 



while at the boundary 

oP, 



dz' 



(66) 






ap.-i , „, ap.-i ap.- 



+ m 



n 



s~2 



dn " dy ' "" dx '" dz 

These are valid for s = 2, 3, 4, ... . We have also 



(67) 



while from (63) we have 



V'Pi 


CVJo 

~ dz' 


3Pi 

tin 


dy 


V'Po 


dva dtiQ 
dx dy 


aPo 


Ivo—tatiQ. 



3Po 

dx 



4- "^Wq 



m 



(660 



(670 



(66^0 
(67'0 



The functions P^ are thus determined successively by the solution of a 
sequence of Neumann's problems. They depend only on the initial motion of 
the liquid and the size and shape of the boundary. 
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We will now express the velocity of the liquid at any time in terms of the 
sequence 1\, 

CO /yS + 2n 

Take ,^^=^^^rrj^,^^ (68) 



Z' 



we have 7s {z) + 7^+2 (^) — ~ , (69) 

|7.(^)=7.-i(2)> (s>0), (70) 

and we notice that 

7q (z) — cos z, 71 {£) = sin z. 

Multiply the first equation in (63) by 7s+i(2ft)^), the second by 7s+2(2ft)i^), 

and add for s = 0, 1, 2,.... The coefficient of ii^ is (2&)^)^/s!, by virtue of 

(69), for s>0, while the coefficient of Wq is 1— 7o(2&)2^). Also, the coefficient 

of '^s is zero for s>0, while the coefficient of ^g^ is —^xi^mf). We therefore 

obtain, on using (61), 

g GO . g coi^ 

'zt = ^to7o(2«0 + %7l(2&>0--5^ 2 Ps7s+i (2(»0--^ S Ps7s+2(2ft)0. (71) 

Similarly, we may obtain 

Q GO 3 '^ 

'^ = '?^o7o(2ft}^)~-^io7i(2wO — ^ 2 l\7s+i(2&)0 + ^ 2 Ps7s+2(2ft)0' (72) 

OJ/ s^O OX s=Q 

10. Methods similar to those of the preceding section may be used to discuss 
the motion of a body which moves irrotationally relative to the rotating axes. 
By this we mean that the body remains . fixed relative to a set of axes which 
themselves always remain parallel to the fundamental rotating axes, but 
whose origin moves in any manner. The most general motion of a sphere 
may be regarded as of this type, when points on its surface are not required 
to be identified. 

These two sets of rotating axes we shall refer to as the fundamental and 
the secondary axes. 

Let u, V, w denote the components at time t of the velocity of the liquid 
relative to the fundamental axes, at a point whose co-ordinates relative to the 
secondary axes are x, y, z. Then, the motion being small, the equations (58) 
remain valid, but we have now 

P = p/p + r2-|a)V, (74) 

m being the distance from the axis of rotation. 

Let V denote the outward normal velocity of the boundary relative to the 
fundamental axes. Then the boundary condition will be 

III -\- mv ~\" niv = Y. 
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Suppose now the expressions (61) and (62) to be valid, at least for a 
certain range of ty together with 

V=iv/-^'. (75) 

5 = 0^- 

Then;^we shall have (63), (65)/ (66), {Q^'), (66'0, as before, but instead of 

(64) we shall have 

lUs-\-mVs-\-niVs = Ys, (76) 

and therefore 

-^-~ — — v^_^^_- V s-i — t — K — -ton—T: 01 —7s — ■} \i i ) 

on dy ox oz 

|i = _V, -Yo-l^ + m ^^ + nwo, (77') 

on oy ox 



'hn 



'Yi + lvQ-miiQ, ' (77") 



for s> 1, instead of (67), (67')? (67''X respectively. 

The formulae (71), (72), (73) remain valid in the present problem. 

We shall now show that the functions P^, for the most general state of 
motion of the type we are considering, provided only that the initial relative 
motion is irrotational, may be derived from those for constant translations 
parallel to each of the fundamental axes, suddenly generated from rest. 

When the secondary axes have a unit translation parallel to Ox of the 
fundamental axes, suddenly generated at time t = 0, and then maintained, 
let P, be denoted by X^. In this case 

Yo=/, V. = Cs>0), (78) 

and the initial motion will be irrotational. Let X-i denote the velocity- 
potential of the initial motion, so that 

V'X-i = 0, ^'= -h (79) 

3X_i 8X-1 3X_i' /oA\ 

no— ~, Vot=—---J:, wo= ^. (80) 

071 ay oz 

Similarly, for unit translations parallel to Oy and Oz, similarly generated 
and maintained, let P^ be denoted by Y^ and Z^ respectively. 

Suppose now that the velocity of the origin of the secondary axes has 
components 

^ (2o)ty ^^(2coty ^ (2o)ty .^.. 

s = ^ ' s = Si 5-0 ^i 

Then it is easily seen from an inspection of the conditions that we shall 
have 

s+l 

x s = Z/ yCtji^s-'n'T' ^n^ s—n'T Cn^s—n)} V^-^/ 

n = 
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for 5= —1,0,1,2,3,..., where P-i will give the velocity-potential of the 
initial relative motion. 

It is now an easy matter to write down the equations of motion of the 
body, if the motion is known to he of the type we are considering. 

As has already been mentioned, it is proposed in a future paper to give an 
account of the small motion of a sphere in a form which is numerically 
intelligible. The expansions used in this particular case have been proved to 
converge for all finite values of t. 
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At the present time there is no generally accepted theory of colour vision, 
and there is considerable divergence of opinion as to the tests which should 
be applied for colour blindness. The Young-Helmholtz theory certainly 
explains the facts of colour mixing, but it is doubtful whether we can deduce 
from these facts the existence of three primary colour sensations. 

For, if we represent colours by particles on a plane, making the masses of 
the particles proportional to the luminosities of the colours, and represent 
mixtures of these colours by the rule for determining the centre of gravity of 
the particles, as is done on Newton's colour diagram, then, if any triangle be 
drawn enclosing all the possible points, its corners may be taken as 
representing the three primary colours, and these corners may take up an 
infinite number of positions. The usual procedure is to make two sides of 
this triangle fit the curve representing the spectrum as closely as possible, and 
take two corners of the triangle on the ends of the curve. Two of the 
primary colours are then definite colours of the spectrum, while the third 
is a sensation which can never be produced pure. But there is a certain 
amount of arbitrariness about this proceeding. 

Also there is no doubt that the belief in primary colours is partly due to 
the experiments in mixing pigments which we are all familiar with from 
childhood, and rests rather on an inveterate tradition than on reason. Just 
as there were four elements, fire, water, earth and air, so there ought to be 
primary colours. I have several times demonstrated to schoolboys, who 



